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ENERGY EQUATION APPROXIMATIONS IN FEBID MECHANICS 
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By Arthur W. Goldstein 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

There is some confusion in the literature of fluid mechanics-*-* ^ in re- 
gard to the correct form of the energy equation for the study of the flow 
of nearly incompressible fluids . Several forms of the energy equation and 
their use are therefore discussed in this note. 

The general energy equation for a fluid, when the chemical composition 
is fixed, may be expressed in terms of the entropy S, the enthalpy H, or 
the internal energy U as 

pTS = pH - p = pU + ppv = J (l) 


where 

p = l/v, density 
T = temperature 
p = pressure 

J = net rate of heat addition per unit volume, including conduction ( -V ' kVT ) , 
dissipation ($), heat sources (pQ) 

The dot indicates the substantial time derivative, e.g., 

U = dU/dt ■ dU/dt + V * VTJ 


When the energy equation is used in this form, no question arises, but when 
H or U is expressed in terms of T, an error can appear by neglecting the 
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The commonly used coefficients needed for the discussion are: 
Cp = (SH/dT)p = H{dS/dT /p, specific heat, constant pressure 
C y = (3U/dT) v = T(dS/dT) v , specific heat, constant volume 
a = (dv/3^)p/v, coefficient of thermal expansion 
Pj = -(3v/3p)m/v, coefficient of compressibility 
P = (3p/3T) y /p, coefficient of tension 

’The existence of an equation of state between p, v, T implies 

a ** PP^p (Eq. (l.05)) 3 
The reciprocity relation is 

(stl - (If), 


0p 


(Eq. (4.23)) 3 


or 


(H) 2 - -($$ v - 


-av (Eq. (4.46))^ 


Also 


( 2 ) 


(3) 


(4) 


Oy = aBpvT 

(Eq. (4.27)) 3 

(5) 

v), then 



H = G p T + 

(l - aT)vp 

(6) 

O = GyT “ ' 

(l - PT)pv 

(7) 


and with the above definitions, the energy equation may be stated in the two 
forms 

<J » pGyT + pTppv (8a) 

«T => pCpT - aTp (8b) 

For a perfect gas, aT = pi = 1, yielding the familiar forms. 


E - G p T, U - GyT, J = pC p T - p = pC v T + ppv 


(9) 
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With liquids it would seem reasonable to neglect the terms involving v 
in Eqs. (7) and (8a). This, however, is not generally permissible because 
0 can become very large . For liquids of very small compressibility ( 0p 
small), the value of p( = a/Ppp) may be large with a consequence that the 
term 0Tppv in (8a) and 0Tpv in (7) might not be negligibly small. The 
kind of magnitudes involved may be judged from the data of the table: 


Liquid 


V, 

a. 

01* 



3, 

1 

ft 

o 


°C 

atm 

1 

1 

cm^ 

Joule 

l 

C-n 




OC 

atm 

gm 

gm-°C 

°C 

P 

HgO 

220 

50 

0.001546 

102.8X10" 6 

1.187 

4.26 

0.301 

0.324 

Hg 

100 

1.0 

0.000180 

4.26X10 -6 

0.075 

0.137 

42.2 

0.157 


The table would suggest that for many liquids 0T » 1. Eq. (7) may 
then be expressed as U « C V T + 0Tpv = pTS = J to the effect that all heat 
addition increases internal energy and none is available for external work. 
The energy equations(8a,b) do not permit such general simplification inde- 
pendently of the particular problem* For example, in a problem where p/p 
is of the same order as T/T, the ratio of the term aTp to pCLT in 
Eq. (8b) is of the order of ap/pCp, which for water and mercury at condi- 
tions cited in the table has values of 2X10“^ and 10"^, respectively. For 
these fluids under such circumstances, Eq (8b) becomes 

J ~ pC T (10) 

0n the other hand, for^ isentropic flows (such as an accoustic wave), J is 
small relative to both pCpT and aTp, and both these terms must be retained. 

In many flow problems, the use of Eq. (8b) is preferable to (8a) because 
the ranges of p and T are indicated by the boundary conditions, ma>l n g 
it possible to decide in advance when either term may be omitted. 
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It is often assumed for liquids that Cp «= (e.g., Ref. 3, p. 4l) 

because of the factor av in Eq. ( 5 ). The effect of the additional factor 
3pT may render this assumption inadequate in many cases as shown in the 
table. In cases where (8b) may be simplified to (lG), the corresponding 
simplification of (8a.), namely J = C v p'i', is not admissible unless the 
liquid satisfies the additional condition Cp « C y , which cannot be assumed 
for all condensed fluids. The particular fluid and problem together will 
determine whether (8a) or (8b) yields the simplest equation. 

The approximation used by Millsaps and Pohlhausen 1 is J = pC v T for 
water at 30° C and unstated pressure, which will be assumed to be 1 atm. 

For these conditions, a = 0.0003/°C, (3p = 48xl0" 6 /atm, from which 
(Cp - CyJ/Cp = 0.014. Since the effect of p is small, the correct 
approximation is J = pCpT. The discrepancy in the energy equation is 
minor because of the small value of (C p - C v )/Cp. The same equation is 
also used for argon where the error is not small, since (C p - C v )/C p = 0.4. 
The proper approximation with small p is J = pCpi as was noted by 
Qstrach and Thornton. ^ 

Hidalgo 2 uses the approximation J = pCpT - p for liquid glass. In 
principle, the equation is not correct, since, if the p term is to be 
retained, it should be retained as aTp (eq. (8b)). In the problem investi- 
gated, p turns out to be negligibly small, so that a physically significant 
error would develop only if a,T » 1. For glasses in the temperature range 
of the problem, aT is not large, so that no substantial numerical errors 


result from this cause 
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